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Due to a need for improved reliability and performance in aerospace systems, there is 
increased interest in the use of adaptive control or other nonlinear, time-varying control 
designs in aerospace vehicles. While such techniques are built on Lyapunov stability theory, 
they lack an accompanying set of metrics for the assessment of stability margins such as the 
classical gain and phase margins used in linear time-invariant systems. Such metrics must 
both be physically meaningful and permit the user to draw conclusions in a straightforward 
fashion. 

We present in this paper a roadmap to the development of metrics appropriate to 
nonlinear, time-varying systems. We also present two case studies in which frozen-time 
gain and phase margins incorrectly predict stability or instability. We then present a 
multi-resolution analysis approach that permits on-line real-time stability assessment of 
nonlinear systems. 


I. Introduction 

We present in this paper the results of our preliminary investigation into the development and use of 
stability metrics for stability analysis of adaptive control in aerospace vehicles. The primary outcome of this 
effort is the development of a road map toward the development and implementation of practical stability 
metrics. Included in this paper are: an overview of the long-term goals, objectives, and expected outcomes 
related to the development of practical stability metrics; §11; and a description of a potential technical 
approach for the development of stability metrics, §111. This description begins with a literature review 
of existing metrics, followed by a description of technology gaps and approaches to fill those gaps. We 
then present preliminary results in our investigation that indicate the failure of current stability metrics 
(gain, phase margins). These same results also indicate the potential for a multi-resolution analysis (wavelet 
decomposition). We conclude this discussion with a formal road map toward the development of stability 
metrics and associated software implementation. 

II. Project overview 

The effort detailed in this paper is the first step to a much larger, long term effort.. The goal of this long- 
term effort is to develop practical theoretic metrics for nonlinear, time-varying systems to facilitate stability 
and performance analysis and design methods for adaptive/intelligent reconfigurable flight control. These 
metrics must meet several objectives, which we enumerate in §A. The outcomes of our road map for the 
development of these metrics are enumerated in §B. 

A. Objectives 

Theoretical stability metrics within this project should achieve the following objectives. 
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1. Applicable to complex nonlinear, time- varying systems with arbitrary order such as a flight control 
system; 

2. Capable of addressing time-scale separation inherent in different actuators, as well as the nonlinear 
dynamic nature of the actuator response; 

3. Capable of addressing stability and performance robustness as well as transient performance relating 
the metrics to classical measures such as gain margin, phase margin; 

4. Capable of facilitating fault prognosis and diagnosis; 

5. Capable of facilitating multi-objective direct, indirect, and hybrid (direct and indirect) adaptive control 
methods to mitigate multiple faults, failures, and damage conditions under uncertain (and potentially 
deteriorating) conditions. 

6. Appropriate for use in flight software verification and validation (V&V). 

7. Capable of balancing speed of response with risk reduction while taking into account the human 
operator. 

B. Expected Outcomes 

Our long-term road map in pursuit of the above objectives is expected to yield the following outcomes: 

1. A necessity stability metric called Departure Index (DI) based on Lyapunov’s first method; 

2. A sufficiency stability metric for equilibrium stability called Stability Margin Index (SMI) based on 
Lyapunov’s second method; 

3. A stability fault warning indicator called Stability Fault Warning Indicator (SFWI) based on Lya- 
punov’s second method; 

4. A sufficiency stability metric called Passivity Margin (PaM) that is the nonlinear, time-varying coun- 
terpart to the linear time-invariant phase margin (PhM) based on passivity theory; 

5. A sufficiency stability metric called ^ 2 -Gain Margin (^ 2 -GM) that is the nonlinear, time-varying 
counterpart to the linear time-invariant gain margin (GM) based on the small gain theorem and 
passivity theory. 

6. Software implementation and simulation demonstration of the use and effectiveness of the proposed 
stability metrics in controller adaptation and reconfiguration (A&R), fault prediction and detection 
(P&D) and flight software V&V. 

III. Technical background 

The development of intelligent and adaptive flight control software for resilient aerospace vehicles is a 
task that is both challenging and essential. This conclusion is based upon several observations: 

• Vehicle models are inherently nonlinear, time- varying and complex, therefore model simplification is 
often used in controller design. 

• The vehicle model, environment, and condition, (component conditions may be nominal/degraded/faulty) 
are uncertain. 

• Vehicle tolerances to off-nominal operation (e.g., stall, side slip, spin) may be limited. 

• Vehicle control authority is limited and may degrade in off-nominal flight conditions. 
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As a result, control flight software must be capable of responding to changing and/or degrading vehicle con- 
ditions while taking into account the uncertainties in the vehicle dynamic model and the flight environment. 
The design of control algorithms that achieve these goals fall into two general categories: robust control and 
adaptive control. 

Robust control refers to a control system's ability to maintain performance and/or stability in the presence 
of modeling errors, whether parametric, dynamic, or both. In the linear time-invariant setting, single- 
input single output systems robustness is commonly assessed in terms of loop gain frequency response 
characteristics, based on Nyquist’s theorem. The simplest and most common quantitative measures in 
this setting are gain and phase margins. More recently, state-space treatment of robustness includes such 
theoretical constructs as 'H- x norms and optimal control, 1 p- analysis and synthesis. 2 which provide measures 
of robustness with respect to specified uncertainty models and connections. Numerous metrics have been 
defined in an attempt to identify robustness margins, such as the graph metric, 3 the gap metric, 4 the point- 
wise gap metric, 5 and the u — gap metric. 6 many of which have been extended to nonlinear systems. These 
results tend to be conservative, similar to results based on the small gain theorem. 7 

Robust control techniques, particularly those based on linear time-invariant (LTI) system theory, tradeoff 
system performance for robustness in a deterministic, worst-case optimal manner. They often treat known 
nonlinear, time-varying (NLTV) dynamics as uncertainties, therefore the resulting control systems, when 
applied to the NLTV plant, are often either marginal or overly conservative in performance and stability 
robustness. On the other hand, adaptive control approaches seek optimality in both performance and 
robustness in response to system behavior change in real time. In principle the adaptive control approaches 
should achieve improved performance over the LTI robust control methods. However, Cusumano and Poolla 8 
discuss capabilities and limitations of adaptation (NLTV control) vs robust (LTI) control in terms of the 
nature of system uncertainty: parametric uncertainty yields advantages to NLTV control, but dynamic 
uncertainty reduces that advantage so that, under some conditions, the achievable performance of robust 
control will approach that of NLTV control. One cause for this limitation is that most adaptation laws 
are synthesized based on LTI stability criteria, which do not apply in general to NLTV systems. Many 
nonlinear adaptive control techniques rely on feedback linearization to cancel the nonlinearity in order to 
use the LTI stability criteria. However, the cancellation is model based and via feedback, therefore it is 
sensitive to dynamic uncertainty (singular perturbation). Another challenge faced by adaptive control is the 
difficulty in software V&V for nondeterministic behaviors due to on-line adaptation. Some adaptive and fault 
tolerant control algorithms also require vehicle failure information form parameter changes to remodeling of 
the vehicle. 

As noted above, most of the LTI-theory based gain- scheduling, robust and adaptive control techniques for 
NLTV systems rely on a frozen-time extension of the LTI stability criteria, either based on the eigenvalues 
(poles) or the gain and phase margins of the linearized system at a fixed time. Although counter examples 
have been provided to show that the stability criterion based on the negative definiteness of the frozen- 
time eigenvalues is neither sufficient 9 nor necessary 10 for the stability of LTV systems a , the success over 
almost half a century of such control techniques in many nontrivial industrial applications, such as flight 
control, seems to suggest that these counter-examples may be merely of academic interest. In fact, the 
widespread success of frozen-time based control techniques can largely be attributed to a (rigorously proven) 
slowly varying constraint, under which the frozen-time eigenvalue stability provides a sufficiency criterion for 
stability. As the demands grow for agility and fast transition through wide operating conditions, the frozen- 
time based control techniques face increasing challenges. Indeed, the modern robust control techniques were 
developed in the hope of expanding LTI based techniques in the presence of NLTV dynamics that are, for 
the convenience of the LTI control techniques, purposely ignored in frozen-time, frozen-state modeling. 

Here we present two practical examples of NLTV systems to show that blind trust and application 
of frozen-time analysis may result in an unreliable stability assessment. Both examples can be modeled as 
pendulums with time-depending damping and stiffness coefficients. Thus, they have broad implications since 
many practical dynamical systems can be decomposed into sets of second-order modes. These examples are 
used later in this paper to demonstrate the effectiveness of our proposed stability metrics based on Lyapunov’s 
first method. 

Example 1. Playground Swing A playground swing can be modeled as a pendulum with a time- varying 
length as the player pumps the swing up by moving the torso up and down, as shown in Figure 1. The 

a See also, 11 where the original Marcus-Yamabe and Wu equations are parameterized to demonstrate the parameter set 
wherein the frozen-time stability criterion fails to be sufficient or necessary for stability assessment. 
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Figure 1. Playground swing: m = 20 (kg), g = 9.81( in/s 2 ), c = 20 (kgm 2 /s), Lq = 3.0 (in), —0.2 < l(t) < 0.2 (m). 


Table 1. Frozen-time analysis of the playground swing (10 cases) 


Frozen Length 

Frozen-time Poles 

Frozen-time Margins 

L+l(t) 

c 

UJ 

Gain Margin 

Phase Margin 

2.800 

1.942 

0.038 

oo 

13.755 

2.844 

1.910 

0.036 

oo 

13.217 

2.889 

1.879 

0.035 

oo 

12.709 

2.933 

1.850 

0.033 

oo 

12.228 

2.978 

1.822 

0.031 

oo 

11.773 

3.022 

1.795 

0.030 

oo 

11.342 

3.067 

1.769 

0.029 

oo 

10.932 

3.111 

1.745 

0.028 

oo 

10.543 

3.156 

1.721 

0.027 

oo 

10.173 

3.200 

1.700 

0.026 

oo 

9.821 


“frozen-time” equation of motion is given by fixing the length of the pendulum at every moment in time as 

c g 


9 + 


m[L + l(t)f 9+ [£+((*)] 


sm(9) = 0 


( 1 ) 


Clearly, for every fixed length /(f), the equilibrium point (0(f), 9(t.)) = (0, 0) is asymptotically stable. Figure 
2 shows the impulse responses and frequency responses of the frozen-time model (1) for ten different fixed 
lengths. Table 1 shows the frozen-time eigenvalues and the frozen-time gain and phase margins. All plots 
and data indicate asymptotic stability. However, the correct equation of motion is given by the principle of 
conservation of angular momentum as 

(i + ;(i) + m[L + i(t)Y) 6+ [z + z(t)] smW = 0 (2) 

Using the parameters shown in Figure 1 and a pumping motion /(f) = 0.2sat(— 9(t)9(t)), the null equilibrium 
(0, 0) can be destabilized with a sizable damping coefficient c > 0, as shown in Figure 3. 
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Figure 2. Frozen-time analysis of the playground swing (10 cases) 


This example convincingly shows that the frozen-time stability criterion fails to be a sufficiency condition 
for stability of LTV and linearizable NLTV systems, and the frozen-time gain margin and phase margin are 
unreliable stability metrics for LTV and linearizable NLTV systems. Noticing that the frozen-time model 
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(a) Playground swing angle 9{t) destabilized by pumping motion 
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(b) Swing length L + l(t) as a result of the pumping motion 


Figure 3. Nonlinear simulation of the playground swing model (2). 


(1) ignored the l(t) term, which becomes negative at times, one may argue that if the frozen-time analysis 
is applied to the correct model (2), the instability may have been correctly predicted. The next practical 
example convincingly shows that, to the contrary, a LTV system may be asymptotically stable even if its 
frozen-time model is unstable at times, thereby showing that the frozen-time stability criterion fail to be a 


6 of 25 


American Institute of Aeronautics and Astronautics 



Figure 4. Jigsaw mounted inverted pendulum: m = 0.5 (kg), g = 9.81 (m/s 2 ), c = 0.01 (kg • m 2 /s), L = 0.5 (m), 
y(t) = Acos(uif), where A = 0.01 (m), = 314 (rad/s) 


necessity condition for LTV and NLTV systems. 

Example 2. (Inverted Pendulum on a Vertically Vibrating Base) An inverted pendulum mounted on a 
vertically vibrating base, such as a jigsaw, is shown in Figure 4. The “frozen-time” equation of motion for 
the inverted pendulum on a vertically vibrating base is no difference from an inverted pendulum mounted 
on a fixed base 


“ + ^“ + (-f) sin(Q) = 0 (3) 

Therefore the upright equilibrium point (a(£), a(t)) = (0, 0)is unstable. However, the correct equation of 
motion is given by the conservation of angular momentum as 

“ + ^“ + (“i + ^ 1 ) sin(Q) = 0 (4) 

where y(t) is the vertical displacement of the base. Now we consider the frozen-time analysis with the effect 
of the moving base with ten different acceleration values y(t) for y(t) = 0.01 cos(314£), which amounts to 
3,000 strokes per minute with a stroke displacement of 0.02 meters. The frozen-time eigenvalues and gain, 
phase margins are shown in Table 2. The corresponding impulse responses and frequency responses are 
shown in Figure 5. These results show that in five out of the ten cases the frozen-time criterion predicts 
instability. Thus it seems reasonable to predict the system would be unstable. 

By letting y(r) = A cos(2t) and d/d r = — d/dt , where ui is the angular frequency of the base vibration, 
and t = tot is the normalized time, equation (4) linearizes to a damped Mathieu equation 

x + 2£ x + (a — 2q cos(2r)) x = 0 , (5) 

whose stability domain is well known 12 and is given by Figure 6. It is seen that under a certain parameter 
combinations of a, q, the upright equilibrium can be stabilized. Figures 7 shows the simulation result for 
the parameters given in Figure 4, which indeed stabilizes the null equilibrium. This example convincingly 
demonstrates the failure of frozen-time criteria as a necessity condition for stability. 

Stability and robustness in nonlinear systems are typically addressed in terms of Lyapunov analysis, 
passivity and the small gain theorem. 1 '* 15 Very roughly speaking, the small gain theorem corresponds 


7 of 25 


American Institute of Aeronautics and Astronautics 


Table 2. Frozen-time analysis of the jigsaw mounted inverted pendulum (10 cases) 


Frozen Base Accel. 

Frozen-time Poles 

Frozen-time Margins 

V (m/s 2 ) 

Ai 

A 2 

Gain Margin 

Phase Margin 

-985.96 

-44.6667 

44.5867 

N/A 

N/A 

-766.86 

-39.4524 

39.3724 

N/A 

N/A 

-547.76 

-33.4336 

33.3536 

N/A 

N/A 

-328.65 

-26.0578 

25.9778 

N/A 

N/A 

-109.55 

-15.4907 

15.4107 

N/A 

N/A 

109.55 

-0.04+ 14.1238; 

-0.04- 14.1238; 

0.9801 

0.0065 

328.65 

—0.04 + 25.2524; 

-0.04 - 25.2524; 

0.9362 

0.0116 

547.76 

-0.04 + 32.8008i 

-0.04 - 32.8008; 

0.8924 

0.015 

766.86 

-0.04 + 38.9114; 

-0.04- 38.9114; 

0.8486 

0.0178 

985.96 

-0.04 + 44.1848; 

-0.04-44.1848; 

0.8048 

0.0203 


to the concept of a gain margin and passivity corresponds to the concept of phase margin. However, it 
appears that these theoretical stability criteria have not been reduced to practical measures that can be 
used in control A&R, fault P&D and control software V&V. Most recently, Zhu et al 16 demonstrated the 
use of Lyapunov’s first, method in the development of a simple metric for analysis and on-line detection of 
instability, where they used an empirical technique to determine nonlinear measures analogous to gain and 
phase (delay) margins. The authors prior work in 16 also alludes to the development of additional practical 
stability metrics using Lyapunov second method, passivity and L-stability theory under a NASA contract 
for the NRA8-30 Space Launch Vehicle Program that was subsequently aborted. 

IV. Approach 

The classical stability indices such as the PhM and GM for LTI systems rely on Laplace or Fourier 
transforms that are not applicable to LTV and NLTV system. Therefore, in general these concepts cannot 
be extended to LTV and NLTV systems using frequency domain or frozen-time approach. Our long term 
goal is to develop practical theoretic stability metrics for NLTV system (including LTV systems as a subset) 
based on Lyapunov’s first and second methods for stability assessment, 17 and the modern nonlinear stability 
theories such as the passivity theory and the ^-stability theory (see, e.g. Khalil 13 ) . 

1. Lyapunov Methods for Assessing Stability of Motion 
Consider the general NLTV system described by state equations 

x = f(t, x), x(t 0 ) = x 0 , t > 0 (6) 

where / is assumed to be piecewise continuous in t and locally Lipschitz in x. A state vector x e (t) £ ]R n is 
called an equilibrium state of (6) if f(t,x e ) = 0. An equilibrium point x e is said to be isolated if there exists 
a ^ > 0 such that for any other equilibrium point x e of (6), ||x — x e \ | > S. In control systems, (6) may model 
the set-point or tracking error dynamics of a closed- loop system, for which the origin is always an equilibrium 
(the null equilibrium) that we wish to attain and maintain under perturbations to the trajectories of the 
system. Thus, it is of interest to determine the stability of the null equilibrium in the sense of Lyapunov. 
It is noted that when the system is time- varying, which is the case for adaptive control systems, uniformity 
of the stability and asymptotic stability must be considered. In order to achieve inherent robustness to 
regular and singular perturbations of the system (structural stability, or robust stability), and bounded- 
input-bounded-output response to disturbances and noises (or £ 0O -stability), (uniform) exponential stability 
should be designed for and maintained, even if only locally (see, e.g. Khalil 13 Chapters 4, 5, 9, 11). 
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(a) Impulse responses 



Frequency (rad/sec) 


(b) Frequency responses 


Figure 5. Frozen-time analysis of the playground swing (10 cases) 


In addition to defining in a rigorous manner the concept of stability of motion, Lyapunov also devised 
two methods for the assessment of stability of nonlinear motion. These two methods are illustrated below 
with the most elementary cases for the null equilibrium. 
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Figure 6. Stability Domain of the Damped Mathieu Equation 12 . 


1. Lyapunov’s first method (indirect method): define the Lyapunov exponent for a solution x(t) = s (i, xo,£o) 

A(xo,to) = limsup ln||s(t ’ Xo ’ to)l1 

t — * OC t 

Then x(t) — * 0 exponentially as t — > oc if X(xo,to) < 0. Thus, the origin is asymptotically stable at 
to if \(xo,to) < 0 for all ||xo|| < 7 , and it is (uniformly) exponentially stable on [To,oo) if X(xo,to) < 

— c < 0, for some c > 0, for all ||xo|| < 7 and for all to € [To, oc). 

2. Lyapunov’s second method (direct method): for an autonomous (time-invariant) nonlinear system ( 6 ), 
i.e. x = f(x ), the origin is asymptotically stable if one can find a positive definite, continuously 
differentiable function V (x(t)) such that V (x(t)) is negative definite (in a neighborhood of the null 
equilibrium). The function V (x) is called a Lyapunov function for ( 6 ) at the null equilibrium. 

It is noted that the second method for NLTV systems may require a time-dependent Lyapunov function 
V(t, s), in which case three time-independent positive definite functions Wi(x), W 2 (x) and W 3 (x) are 
required such that 


W 1 (x) < V{t,x)<W 2 (x) 
dV dV 

V{t,x(t)) = ~^ + -^f(t,x)<-W 3 (x) 

Moreover, for exponential stability, it is required that Wi(x) = kj ||x|| n for some positive constants a and k{, 
i = 1,2,3. The second method is a sufficiency condition in that the inability to find a Lyapunov function 
does not imply the equilibrium is unstable. 

It is also noted that Lyapunov’s first method is often mistaken as the linearization method in the literature 
and many textbooks. In fact, the first method was only applied by Lyapunov to LTI and LTV systems to 
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Figure 7. Inverted pendulum stabilized by the pumping motion at the base 


derive some general results. The linearization method, which is based on the important equivalence condition 
that the null equilibrium of a NLTV system (6) is exponentially stable if and only if the linearized system is 
exponentially stable, is established using the second method. 


2. Lyapunov First Method Based Stability Metrics 

While the Lyapunov exponent for a particular (tracking error) trajectory x (t) = s (t , x - o, to) is readily com- 
putable in real-time by 


Ao = lim sup 

t — *oo 


r l 1 

i(r) 

dr 

Jt o 1 

x(t) 


i 

1 

o 



for some convenient norm and its negativity can be monitored, a difficulty in using the first method in 
assessing stability is that it requires the solution of the nonlinear system for every possible initial condition. 
However, this does not prevent it from being used as an instability assessment for NLTV systems. Another 
difficulty is that the exponent is defined for t — > oc, while any practical flight mission is of finite duration. 
For time-dependent systems, the system trajectory may diverge, also known as departure, momentarily so 
as to violate some operational constraints, even though it may ultimately converge to zero as t — > oc . In 
order to detect such practical instability, a short-time, or windowed, Lyapunov exponent can be defined as 

rf ||x(t)|| 

At (t) = J ‘- T T x(r> dr 

where T is the width of the moving window over which the exponential departure rate is observed. In the 
sequel, we shall call Ao a Long-term Departure Index (LDI) for the trajectory, and Xp a Short-term Departure 
Index (SDI). Being necessity conditions for stability, the mam utility of this method is to facilitate on-line 
warning of onset of loss of stability. Figure 8 shows an implementation of the LDI/SDI described above in a 
tracking control system, where x(t) is the tracking error. It is noted that the magnitude of the DI \\t | can 
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Figure 8. Implementation of the DI. 


be used as a transient performance indicator as it is the exponential decay rate of the (tracking error) state 
variable. 

The LDI and the SDI have been implemented and the effectiveness of the DIs has been studied under 
NASA NRA8-30 for reusable launch vehicles on a high fidelity 6-DOF simulation platform MAVERIC at 
Marshall Space Flight Center. Very encouraging preliminary results were reported Best et al. 16 In this paper, 
we apply this method to aircraft adaptive flight control systems to verify it’s effectiveness and to address some 
of the issues identified in the previous study. In particular, the method has a technical difficulty in coping 
with trajectory zero- crossings of a single state variable. This problem can be solved in theory using Hilbert 
transform. However, the Hilbert transform is an anti-causal process and cannot be realized in real-time for 
nonstationary (time-varying) processes. Our preliminary suggestion for this method for circumventing this 
difficulty is to use short-time, narrow-band (wavelet) approximations to the Hilbert transform, which can be 
implemented in real-time using a bandpass filter bank. The theoretical study will focus on approximation 
error analysis and validation of the concept and methodology. Use of the modal information provided by 
the wavelet analysis for diagnostic and prognostic of flight condition and vehicle health, and for control 
adaptation will also be investigated. 

3. Lyapunov Second Method Based Stability Metrics 

In addition to assessing stability of an equilibrium state, the second method of Lyapunov also offers a 
(conservative) estimate of the Domain of Attraction with the ensuing Lyapunov function V(t, x), wherein 
trajectories are guaranteed to converge to the equilibrium. The DoA is specified by the interior of a closed 
level contour T c = {x : W\ (x) = C > 0} wherein — W 3 (x) < 0. A practical sufficiency stability metric can be 
defined as the distance of the current (tracking error) state trajectory x(t) from the boundary of the DoA 

dist ( x(t ), r c ) = inf || x(t) — z\\ 


using a suitable norm, which will be termed Stability Margin Index (SMI), which can be monitored with full 
state measurement or estimate. 

This method depends on an a priori estimate of Domain of Attraction (DoA) using a Lyapunov function for 
the closed-loop system, which is usually available for controllers designed using Lyapunov’s second method. 
Otherwise an estimate of the DoA may be obtained using an iterative method to construct a Lyapunov 
function, 18 or by reverse-time simulation method. The limitation of this method is that it is a model based 
method. Thus its accuracy and validity rely on the fidelity of the model, and models for adverse flight 
conditions and vehicle failures need to be developed in order to develop the corresponding estimate of DoA. 

As a byproduct of this method, another index to monitor is the negativity of the time derivative of the 
Lyapunov function along the current trajectory. Although the negativity alone is neither a necessary nor a 
sufficient condition for stability, loss of negativity in an a priori validated DoA can be used as an indicator 


for abnormal flight conditions or vehicle health. Moreover, the magnitude 


\V(x)\ 


is an indicator of transient 


performance when V (x) < 0, and the magnitude \V (x) | is an indicator of steady-state performance indicator. 
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4 - Assessment of Input- Output Stability and Stability of Feedback Loops 
Consider a NLTV control system described by 

x(t)=f(x(t),u(t), 0 (t)) ( 7 ) 


y(t) =h(x(t),u(t),e(t)) 

where x(t) £ R n , u(t) £ R ; , y(t) £ R m are the state, input and output vectors , respectively, and 0 (t) £ R p is a 
piecewise continuous time- varying parameter vector representing time-dependent flight condition parameters 
such as Mach and dynamic pressure, system parameters such as aerodynamic coefficients, and adaptive 
controller parameters (gains). The function /is locally Lipschitz and h is continuous. The NLTV system 
( 7 ) defines an input-output (I/O) map y — Hu , which can be used to model the subsystems in a feedback 
loop, such as the plant, actuators, sensors and controllers. It is noted that the overall closed-loop system 
with a reference input is itself an I/O map of the form ( 7 ). 

It is of practical interest to consider whether the output response y = Hu to every norm-bounded input 
u is norm-bounded for some suitable norms. If so, the I/O map H is said to be input-output stable with 
respect to the norms of the I/O signals. Commonly used norms in this context are the extended T p norms 
(for definition, see e.g. Khalil 13 pp. 195 - 197 ). In particular, if for every C p norm-bounded input signal w, the 
system ( 7 ) satisfies 



< 7 


for some positive constant 7 then I/O map H is said to be (zero-bias) finite-gain Testable, and the supremum 
of 7 taken over the entire extended T p -space of the input signals is called the gain of the I/O map. We will 
subsequently use 7 to denote the gain of the I/O map H. Two commonly used C p norms are £2 and T^, 
representing the energy and peak-magnitude of the signal over a time-interval. Thus finite-gain /^-stability 
signifies bounded-energy I/O stability, and the finite-gain Toe-stability is the generalization of the Bounded- 
Input-Bound-Output (BIBO) stability for LTI systems, the latter is equivalent to (exponential) asymptotic 
stability of a LTI system. 

It is noted that, unlike LTI systems, for nonlinear systems, the null equilibrium in the (tracking error) 
state space being merely asymptotically stable does not guarantee I/O stability. The practical implication 
is that a bounded disturbance or noise input may cause the system output departure. An important result 
in the £ p -stability theory is that, under some mild conditions on / and h, exponential stability implies 
(small-signal) finite-gain Testability of the I/O map. 

A11 important Testability result is the small- gain theorem for a feedback loop shown in Figure 9 , which 
states that the I/O map from any input to any output of the feedback configuration is finite-gain T p -stable 
if the loop gain 7 = 7172 < 1 , where 71, 72 are the T p -gains for the subsystems Si, S2, respectively. This 
result can be viewed as a nonlinear extension of the LTI gain margin concept where for stability the frequency 
response loop gain is required to be less than 1 above the crossover frequency. Thus the margin 7 — 1 can 
be viewed as an T p -gain margin. 

There are two practical difficulties in applying the T p -gain margin concept in real-time stability assess- 
ment. First, it is impractical to estimate the T p -gain in real-time for NLTV systems. Second, the small-gain 
theorem is an overly restrictive sufficiency condition that precludes high-gain control for improved steady- 
state precision. In fact, for LTI systems, the loop gain can be greater than unity below the crossover 
frequency, provided that the phase margin is positive. Both of these difficulties can be circumvented by the 
passivity theory for NLTV systems. 

In a nutshell, the passivity theory is an energy based approach. A system is passive if the power dissipation 
by the system is greater than or equal to the rate of exchange of stored energy in the system, which can be 
stated as 

u{t) T y(t) > V(t,x{t)) = + yy/(T^) 

where V (t, x) is a continuously differentiable positive definite function known as a storage function. The 
system is said to be strictly passive if the inequality holds strictly, in which case the margin 


i/j(x) = w(t) T y(t) — V ( t , x(t)) >0, Wx ^ 0 
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Figure 9. Feedback configuration of two I/O maps. 


is known as the excess of passivity. The excess of passivity for the open-loop I/O map e\ — ► in the 

feedback configuration shown in Figure 9 can be viewed as a NLTV extension of the concept of phase margin 
for LTI systems, because the cascade of two strictly passive LTI systems has a positive phase margin. The 
passivity concept can be extended by considering I/O maps that can be rendered passive via a local feedback 

u(t) T y(t) > V (t, x(t)) + y T p(y), y T p(y ) >0, \/y ± 0 


which is called output strictly passive. The main result of the passivity theory is that the closed-loop I/O 
map is ^-stable, and the null equilibrium of the zero-input closed-loop system is asymptotically stable if 
both subsystems are strictly passive, or one or both are output strictly passive with an additional requirement 
of being zero-state observable, a mild condition that can be satisfied with practical sensors. In addition, 
an estimate of the ^ 2 -gain can be obtained from real-time measurement of the stored energy in the state 
variables, and the output signals (see, e.g., [13, Ch. 6]. 


5. Passivity Theory Based Stability Metrics 


We propose to use the excess of passivity of the open-loop I/O map, called Passivity Margin (PaM), as 
a NLTV extension to the LTI phase margin concept, and use the margin of the ^ 2 -loop-gain from unity, 
called £ 2 -Gain Margin (/I 2 -GM) as a NLTV extension to the LTI gain margin concept. Both margins can 
be monitored using real-time measurement or estimate of the state variables, and input and output signals. 
Parasitic energy stored in the neglect fast modes, e.g. actuator dynamics, structural or aeroelasticmodes, will 
be accommodated by the passivity margins being monitored. Additional modes can be added to the storage 
function as needed through experimentation to improve fidelity. A sudden loss of passivity margin would 
imply onset of adverse flight condition, excitation of unmodeled modes or vehicle faults that either inject 
additional energy into the system, or prevent them from being damped out. Thus, the real-time measurement 
of these margins can be used for fault P&D and control A&R. In particular, the energy measurement can 
be done in divided frequency bands to aid fault analysis and adaptive control law synthesis using, for 
example, the time-varying bandwidth filters (TVBF) or time-varying notch filters (TVNF) developed by 


Zhu et. al. 19,20 Similar to the Lyapunov second method, the magnitude I V (x) 


is an indicator of transient 


performance when V ( x ) < 0, and the magnitude \V (x) | is an indicator of steady-state performance indicator. 
Moreover, the integral of w(t) T y(i) from the control input to the plant output over a time interval is an 
indicator of the control energy consumption over that time interval. 


6. Applications of the Proposed Stability Metrics 

It is noted that none of these methods is by itself a necessary and sufficient condition for stability. Therefore, 
synergistic interpretation of these stability metrics need to be studied and applications of these stability 
metrics in controller A&R, fault P&D and flight software V&V need to be demonstrated. 
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Figure 10. A 6-DOF Flight Controller Using TVNF Adaptation for Flex Mode Suppression 
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Figure 11. Multiplicative adaptive control assist architecture. 


There are two basic approaches to control A&R, direct and indirect. The former adjust the controller 
without explicit plant model parameter estimation or fault detection. One of such techniques is the time- 
varying bandwidth filters (TVBF) or time-varying notch filters (TVNF) for adaptive and fault tolerant 
control developed by Zhu et. al. 21,22 The multi-band DI, PaM and £ 2 -GM can be used to trigger control 
adaptation without explicitly identifying the plant parameters or faults to facilitate effective damping of 
excess kinetic energy and bleeding of excess potential energy. Figure 10 shows a 6-DOF flight controller 
utilizing the TVNF adaptation for a flexible hypersonic vehicle model developed by AFRL Air Vehicle 
Directorate 22 based on multi-band tracking error energy measurement, which can be modified to use the 
multi-band DT or PaM with potentially improved performance and robustness. An alternative approach is 
to use adaptive control involving neural networks as a parameterization of the control law where the adaptive 
portion of the controller may be used in a multiplicative (Figure 11) or additive (Figure 12) architecture. 
Adaptation design is performed based on Lyapunov stability criterion 15,23,24 which may be modified to 
take advantage of the on-line estimated stability margins. For indirect adaptive and fault tolerant control 
schemes, parameter estimation or fault diagnosis will need to be explicitly performed. 

The proposed stability metrics can used for fault P&D in many ways. For example, the multi-band DI 
may provide modal information to indicate possible causes for the loss of stability — high frequency may 
indicate actuator failure, excitation of structural or aeroelastic modes, whereas low fr equency may indicate 
failure in control effectors that hampers damping of vehicle kinetic energy. The multi-band PaMs indicate 
energy build up in different frequency band, which also provide hints to the cause and type of vehicle fault 
or flight condition. The correlation between the stability metrics and vehicle failures can be used by indirect 
adaptive or fault tolerant control schemes. 

Finally, the proposed stability metrics may be useful in flight software V&V in several ways. For example, 
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plant 

command 


Figure 12. Additive adaptive control assist architecture. 


once proving effective, adaptive flight software may be required to implement one or more stability metrics, 
and demonstrate that the software respond to these metrics in a predictable manner. The metrics can also 
be implemented in high fidelity, Monte Carlo simulations to validate adaptive flight software by exhaustive 
simulation test. The stability metrics can be used to specify requirement, or to demonstrate improvement 
in robustness and reliability. 


V. Real-time stability assessment by multi-resolution analysis 


Consider the n-dimensional NLTV tracking error or operating error dynamics 

x = f(t, x ), x(t 0 ) = x 0 , t > t 0 (8) 

where / is continuous in t and Lipschitz in x, and x(t) = 0 is the desired equilibrium state that needs to be 
stabilized. The time-dependence of / may be due to either time-varying parameters or time-varying nominal 
tracking trajectory of the NLTV system. For simplicity, here we assume that (8) can be linearized at the 
operating (null equilibrium) point to a LTV system of the form 


x = A c (t)x , x(t o) = xo, t > to 

where A tc (t) is continuous and bounded, and is in the companion canonical form 


( 9 ) 


0 


A c (t) = 


0 


In - 1 


-ai(t) -a 2 (t) 




where I n - 1 is the (n — 1) x (n — 1) identity matrix. By letting y = x i, the state equation (9) is equivalent 
to an nth-order scalar LTV differential equation 


2/ (n) + (*n{t)y (n l) + ■ ■ ■ + a 2 {t)y + ai(t)y = t) (10) 

y ( - k ~ 1) {to) = xo k , k=l,...n t>t 0 

The general solution to LTV equation (10) can be written as 

y(t)= ^2 C k yk(t)= J2 C ke ft o pk{T)dT ( 11 ) 

fc=l k = 1 

where {yk{t) = exp (f pk(t) dt )}” _ 1 is a set of linearly independent solutions to (10), and {pfcfy)}/! =1 is called 
a parallel differential (PD) spectrum for (10) in which each pk(t) is called a PD -eigenvalue 25 ' 26 for A c (t). 
The PD-eigen values are solutions to an (n — l)th-order NLTV differential equation, called PD- characteristic 
equation. In particular, for n = 2, the PD-character equation is given by the Ricatti equation 13 

p + p 2 + a 2 (t)p + ai(t) =0 (12) 

b The general nth-order PD- characteristic equation is omitted here for brevity, for details see Ref. 25 ’ 26 . 
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Note that for real- valued coefficients Qfc(t), the PD-eigenvalues pk{t) may be complex- valued in order to 
avoid singular solutions. Moreover, the synthesis formula for afc(t) from a PD-spectrum {pk(t)}k—i is given 
by 


a l(*) — Pl(t)p2(t) + 


pl(t)p2(t) - p 2 (t)p 1 (t) 
P2(t)-pi(t) 


= -Pl(t)-p 2 {t) 


P2(t) ~ Pli*) 

P2(t) ~ Pl(t) 


(13) 


The PD-eigenvalue concept is a natural extension of the eigenvalue concept for LTI systems in that, 
when the coefficients a*,(t) = ak ^constant, the PD-character (differential) equation (12) possesses constant 
solutions pk(t) = pk ^constant that are conventional characteristic roots for (10), and eigenvalues for 
A c (t.) = A c , for which (12) reduces to the conventional algebraic characteristic equation, and the synthesis 
formula (13) reduces to the conventional polynomial coefficient synthesis formula. Moreover, stability of the 
LTV system (9) or (10), and consequently the null equilibrium of the NLTV system (8), can be determined 
by the negativity of the real-part of the PD-eigenvalues. In particular, exponential stability is implied by 
R e(pk(t)) < —c < 0, for some constant c > 0, and for all k. 

It is also noted that for each k, Pk(t) can be written as 


Pk{t) = 


■i/M 

Vk{t) 


(14) 


where yk(t) is complex- valued if pk{t) is. Thus pk(t ) has a frequency dimension t -1 , which may be termed 
instantaneous frequency of y(t ). The notion of instantaneous frequency was original introduced for frequency 
modulated (FM) signals. It has been generalized to a much richer concept and extensively studied in the area 
of signal processing for signals generated by time-dependent processes, such as audio and visual signals 27-29 
However, the powerful tool of time-dependent frequency analysis, which is very pertinent to LTV and NLTV 
systems, has not been widely utilized in the area of control systems. 

Figures 13 and 14 are taken from Ref. 28 to illustrate the method of time-dependent frequency analysis 
for single-mode and multi-mode signals, respectively. The plots, called spectrograms, describe the harmonic 
distribution in the signal as function of time. These signals can be produced, respectively, by second-order 
and higher-order LTV systems, and there is an intrinsic relationship between the spectrogram and the PD- 
eigenvalues of the LTV system that can be utilized for real-time stability assessment, as well as LTV system 
identification. 

Consider a complex-valued single-mode signal 


x(t) = Ce jf ^ t)dt 


This signal can be produced by a 2nd-order LTV system 

x - + u 2 (t)x = 0 

uj{t) 


for which pi t2 (t) — ± 7 *u>(£) are PD-eigenvalues that constitutes a PD-spectrum. We now demonstrate the 
relationship between the spectrogram and the PD-eigenvalue by way of an example. Let u j(t) = — 6sin(0.3t), 
and denote by Pa(^) the unit intensity pulse function centered at t = 0 with magnitude A -1 and width A, i.e. 


Pa(*) = 



*e(-#,f] 

elsewhere 


Let x(t) be a function composed of complex sinusoidal splines 

OO 

x(t) = ^ p A (t- kA) e JWfct A 

k= — oc 
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(a) Time trajectory of 9(t) 




loglO(freq) 


20 

time 


(b) Spectrum of 9{t) over [0, 100] 


(c) Spectrogram of 6(t) 


Part (a) shows the time trajectory 9(t) of a playground swing discussed in this paper. The oscillation not only 
diverges in magnitude, but also decrease in frequency. The Fourier spectrum shown in Part (b) does not capture 
either of these important features of the oscillatory mode. But these features are captured by the Spectrogram 
shown in Part (c), see also Figures 3 and 17 and the discussions thereof. 

Figure 13. Spectrogram of a single-modal signal. 



(a) Time trajectory of u(t) 



(b) Spectrum of u(t) over [0, 100] 


25 -v 



loglO(freq) ' 2 0 timB 

(c) Spectrogram of u(t) 


Part (a) shows the tracking error of the body frame velocity u(t) along the ar-axis of a 6DOF vertical take-off and 
landing (VTOL) UAV when it flies a semicircular trajectory. The doublet tracking error is due to the acceleration 
and deceleration at the end points of the semicircular trajectory. The simulation consists of 66th-order dynamics. 
Part (b) shows the Fourier transform of the tracking error signal over the 100-second time interval. The spectrogram 
in Part (c) not only shows the energy concentration in the frequency domain with much higher resolution, but also 
the time when the energy built-up happened, and the exponential convergence of the tracking error. 


Figure 14. Spectrogram of a multi-modal signal. 
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where uik — w(k A). For sufficiently small A, the piecewise constant function 


u}(t) = U>fc, 


t G 




2k +1 
2 


A 


is a piecewise constant approximation of u>(t ), and x(t) is a sinusoidal spline approximation of x{t). Now tak- 
ing the short time Fourier transform (STFT) of x{t) with the kth. analysis (observation) window P&(t — A: A), 
we have 


Xk{ju) = sine 




Figure 15 shows the plots of u/(t), u)(£) and Xk(ju>) for different analysis window sizes. It is noted that the 


collection of 


X k (juj)\ constitute a spectrogram for x(t ), where 


with main lobe bandwidth 


max 


X k {ju) 


= A = 


Xk(ju>k) 


Bm( A) = 4tt/A 


Thus, we can estimate the (imaginary part of) PD-eigenvalue from the modal peaks of the spectrogram at 
t. = A:A. In order to improve the accuracy of the estimate t o(t), it is desirable to use smaller time window size 
A. However, as A decreases, the spectral peak is flattened, making an accurate reading the center frequency 
more difficult. In fact, as A — * 0, 


max 

u > 



0 


Bm{ A) = 4tt/A -► oo 


i.e. the STFT spectrum \Xk(ju)\ approaches the constant zero “centered" at cjfc, as can be seen in Figures 
15(b) and 15(c). However, as A — > 0, Pa (A) — ^ <5(t), it follows 

ui(t) — > uj(t) 


/ +oo 

S(t-r)e if ^ T)dT dr = x(t) 

-oo 

This shows that u>(t) is in fact the precise measurement of the instantaneous frequency at the instant t 
(see Figure 15c). Thus it is highly desirable to develop methods that directly estimate the instantaneous 
frequency from the signal x(£), and then construct an accurate spectrogram from the PD-spectrum. 

The signal x(t) in this example does not decay with time. For signals whose amplitude envelope decays 
(or grows) with time, the PD-eigenvalue will have a real part, which can be estimated by looking at the peak 
amplitude of the spectrogram, as demonstrated in the sequel by the swing and jigsaw pendulum examples 
using the MRA method. 

The imaginary part u>d(t) of the PD-eigen values p(t) = a(t) ± ju!d(t ) can be identified by connecting 
the peak magnitudes of the frequency spectrum at each time t, and the real part cr(t) can be identified by 
calculating the slope of ln(X^,(u;, t)) at each time t. 
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-IfT' 0 


(a) Large Observation Window, A = 2.0, (b) Small Observation Window A = 0.4 

u(t) = — 6sin(0.6 k), 2k - 1 < t < 2k + 1 , Q(t) = -6sin(0.12fc), 0.4fc - 0.2 < t < OAk + 0.2 

= 2sinc(u; -w*), B m ( A) = 2n |Xfc(jw)| = 0.4sinc (0.2(o; - Vk)) , B m ( A) = 107r 



(c) Observation Window A — > 0, kA — * t 
uj(t) = — 6sin(0.3A:A) —* — 6sin(0.3f) = u>(t) 
Xk(jA ) I = Asinc (f (w - u>k)) — ► 0, B m ( A) = f -*oo 


Figure 15. Relationship between the STFT spectrogram and PD-eigenvalues. u>(t) = — 6sin(0.3t). 
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Mother wavelet for stability analysis (a — 5) 



Fourier transform of mother wavelet, a — 5 



Figure 16. Mother wavelet used for stability analysis. 


VI. Case study 


In this section we present a preliminary case study of real-time stability assessment for NLTV systems 
by applying one of the proposed approaches, namely, Lyapunov first method, to the playground swing and 
jigsaw mounted inverted pendulum examples to demonstrate the effectiveness of the proposed approaches 
over the conventional stability metrics. The study is facilitated by multi-resolution (wavelet) analysis (MRA) 
of the tracking or operating error state trajectory, which is a useful technique that is applicable to both 
trajectory based (Lyapunov first) and energy based (Lyapunov second and passivity) stability assessment. 
We demonstrate the utility of this method with the systems examined in Example 1 (pumped swing) and 
Example 2 (jigsaw). 

For this case study, define the family of wavelets 


V>(a, t) 


aei at t € [0, 

0 otherwise 


(15) 


with corresponding Fourier transform 


^(a, uj) =fty{a,t)\ = e^^sinc ^ 


(see Figure 16). Given a signal x(t) we define its wavelet transform 


X^(uj,to) = (iJj(w,t-to),x(t)) = f - t 0 )dt 

J — oo 


(Recall that an inner product in the time domain corresponds to a multiplication operation in the frequency 
domain, with appropriate scaling for phase information.) 

Recall for these systems in Examples 1 and 2 that frozen time stability analysis of the swing model predicts 
that it is stable. However, as seen in Figure 17, this prediction is incorrect. The system response, shown in 
the upper left plot, shows an increasingly large oscillation throughout the first 50 seconds of operation. The 
upper right plot is a multiresolution analysis of the system response, a plot of log 10 \Xjp | as a function of time 
and frequency. The “height" of the peak of the MRA is shown in the lower left plot, and provides a local 
measure of exponential/damping behavior of the system - that is, an assessment of stability. The frequency 
associated with this peak is marked in white in the lower right plot. That is, wavelet analysis permits the 
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Figure 17. Multiresolution analysis analysis of swing behavior. LTI analysis predicts stability; wavelet analysis 
quickly detects instability. 


identification of the time- varying frequency and exponential envelope that describe the instability. Subject to 
the limits of Heisenberg’s uncertainty principle, 30 this MRA identifies time-varying frequency and damping 
characteristics of the system response. 

Conversely, frozen time stability analysis of the jigsaw model Example 2 predicts that it will be unstable. 
However, as seen in the upper left plot in Figure 18, the system response in this simulation appears to behave 
as a damped sinusoid. Multi-resolution analysis (see the upper right plot of Figure 18 reveals the frequency 
ujq of oscillation and, in turn, the slope of log \X^(uio, t)\ reveals that the system can be well approximated 
by a pair of complex conjugate poles with negative real parts. Multi-resolution analysis also detects a nigh 
frequency oscillation in the first few seconds of simulation, which is also visible on the time-domain plot 
(upper left). 

Remark 1. The MRA spectrogram method for real-time stability assessment presented herein can be ex- 
tended to higher order (multi-mode) systems by first isolate the modes using bandpass filters, followed with 
single-mode analysis. This method is applicable to LTV and linearizable NLTV systems. However, for NLTV 
systems, the nonlinearity may cause spectral distortion by creating higher harmonics or sub-harmonics that 
may be mistaken as additional modes. This may cause problem for modal analysis, but it should not be a 
problem for energy based method, as passive nonlinearity does not inject energy into the signal. 

Remark, 2. It is noted that while the spectrogram analysis can be used to assess stability in real-time using 
the PD-eigenvalues estimated from the real-tune measurement of a particular state trajectory x(t), it is not 
reliable to ascertain stability of the null equilibrium based just on this one measurement. There may be 
dynamic modes that are not excited in a particular trajectory. Thus, to assess stability of the equilibrium, 
all modes need to be persistently excited. Moreover, for time-varying systems, it is important to determine 
uniform stability. However, this would require to test the system for infinitely large initial time to. which is 
not practical. But this should not be a problem for engineering applications as all real world systems have 
a finite mission duration, and we may assume the system parameters freeze at the mission terminal time. 

Remark 3. With persistent excitation, the spectrogram analysis can be used to identify a model for a LTV 
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Figure 18. Wavelet analysis of jigsaw behavior. LTI analysis predicts instability; wavelet analysis detects 
stability. 


system using the estimated PD-spectrum. It is imperative to use the correct synthesis formulation to obtain 
the time- varying parameters for the identified LTV system. In fact, any frequency domain LTI spectral identi- 
fication algorithm can be used to estimate the PD-eigen values at frozen time moment from the spectrogram. 
The fundamental difference between the (correct) PD-spectral synthesis and the (unreliable) frozen-time 
method is that the latter uses the algebraic polynomial synthesis formula on time-varying eigenvalues, which 
introduces errors that are dependent on the rate of change of the PD-eigenvalues. This may explain why 
the froze-time based gain-scheduling, robust and adaptive control techniques have been successful for slowly 
varying systems, but are challenged by fast varying systems. 

Remark 4. In the example illustrating the relationship between the spectrogram and PD-eigenvalues, it is 
noted that a large time window (step) yields poor time resolution, but better frequency resolution, while 
a smaller time window provides a better time resolution but poor frequency resolution. This phenomenon 
is known as the uncertainty principle for the STFT. While the MRA method improves the STFT by using 
smaller time-windows for higher frequency harmonics, the ultimate solution, as shown in Figure 15, is 
to directly compute the PD-eigenvalues using equation (14). However, the difficulty in this method is 
that, for complex modes, the imaginary part of the real- valued signal must be estimated. This may be 
accomplished using the Hilbert transform. However, Hilbert transform is an anti-causal operation that 
cannot be applied to real-time signals. Therefore, approximations must be made, a subject that deserves 
future study. Alternatively, a quadrature signal, i.e. another solution trajectory that is linearly independent 
of the measured trajectory can be constructed to facilitate the estimate of PE-eigenvalues using (14). 

Remark 5. The MRA spectrogram real-time stability assessment method presented herein is based on Lya- 
punov’s first method (spectral method). However, the spectrogram can also be used to facilitate stability 
assessment using energy-based stability criteria, such as Lyapunovs second method and passivity methods. 
This is based on the Parseval principle that the energy contained in a signal is the same in its time-domain 
and frequency domain representations. 


23 of 25 


American Institute of Aeronautics and Astronautics 


VII. Conclusions 


The need in aerospace systems for reliable operation in high-performance operating conditions has led 
to increased interest in the use of NLTV (e.g., adaptive or reconfigmable) control laws for these systems. 
Adaptive control laws, while based in Lyapunov stability theory, do not easily fall into existing theory for 
verification and validation (V&V), and thus can be considered a risk in an overall design. The weaknesses 
of classical LTI stability assessment are shown in two simulation examples, one of a playground swing and 
the second of a jigsaw. 

In this paper, we provide a roadmap to the development of stability metrics appropriate for use in V&V 
and in on-board stability assessment. The conceptual basis for these metrics is drawn from well established 
stability concepts such as the small gain theorem, Lyapunov’s first and second methods, and passivity theory. 
We also present a multi-resolution analysis approach that provides real-time assessment of stability in each of 
the example systems. The results of this analysis are closely tied to LTV system concepts (PD-eigen values, 
etc.) and are easily interpreted. 
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